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Abstract 

In the framework of the simplest little Higgs model (SLHM), we perform a comprehensive study 
for the pair productions of the pseudoscalar boson rj and SM-like Higgs boson h at LHC, namely 
gg{bb) — > nn, gg(qq) — > nh and gg(bb) — > hh. These production processes provide a way to probe 
the couplings between Higgs bosons. We find that the cross section of gg —> nn always dominates 
over that of bb —> nn. When the Higgs boson h which mediates these two processes is on-shell, 
their cross sections can reach several thousand fb and several hundred fb, respectively. When the 
intermediate state h is off-shell, those two cross sections are reduced by two orders of magnitude, 
respectively. The cross sections of gg —> nh and qq — > nh are about in the same order of magnitude, 
which can reach O(10 2 fb) for a light n boson. Besides, compared with the SM prediction, the cross 
section of a pair of SM-like Higgs bosons production at LHC can be enhanced sizably. Finally, we 
briefly discuss the observable signatures of nn, nh and hh at the LHC, respectively. 

PACS numbers: 12.60.-i,12.60.Fr,14.80.Ec 
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I. INTRODUCTION 



in 



Little Higgs theory |l[ has been proposed as an interesting solution to the hierarchy prob- 
em. So far various realizations of the little Higgs symmetry structure have been proposed 
-5|, which can be categorized generally into two classes [6|. One class use the product 
group, represented by the littlest Higgs model in which the SM SU(2)l gauge group is 
from the diagonal breaking of two (or more) gauge groups. The other class use the simple 
group, represented by the simplest little Higgs model (SLHM) (4], in which a single larger 
gauge group is broken down to the SM SU{2) L . 

Since these little Higgs models mainly alter the properties of the Higgs boson, hints of 
these models may be unraveled from various Higgs boson processes. The phenomenology of 
Higgs boson in these little Higgs models has been widely studied (t- 14] . In addition to the 
SM-like Higgs boson h, the SLHM predicts a pseudoscalar boson 77, whose mass can be as 
low as O(10 GeV). The constraint from the non-observation in the decay T — > 777 excludes 77 



with mass below 5-7 GeV 



15] . In this paper, we will focus on the pair productions of neutral 



Higgs bosons at LHC in the SLHM, namely gg(bb) — > gg(qq) — > r]h and gg{bb) — > hh. 
These production processes at LHC are very important because they will provide a way 
to probe the couplings of hrji] and hhh, and shed light on the Higgs potential. Further, 
Higg-pair production at LHC may be sensitive to new physics, and it has been studied in 



many new physics model, such as little Higgs models [13|, ]14J, supersymmetric models 
models of universal extra dimensions and left-right twin Higgs model jl8 |. Note that 
gg(bb) — > hh process has been studied in 13|. In order to perform a comprehensive study 
for the pair productions of neutral Higgs bosons at LHC, we reconsider them here. 

This work is organized as follows. In Sec. II we recapitulate the SLHM. In Sec. Ill we 
study ggibb) — > 7777, gg(qq) — > rjh and ggipb) —> hh production processes at LHC, respectively. 
Finally, we give our conclusion in Sec. IV. 



II. SIMPLEST LITTLE HIGGS MODEL 

The SLHM is based on [SU (3) x U (l)x] 2 global symmetry. The gauge symmetry SU (3) x 
U(l)x is broken down to the SM electroweak gauge group by two copies of scalar fields $1 
and $2, which are triplets under the SU(3) with aligned VEVs fi and f 2 . The uneaten five 
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pseudo-Goldstone bosons can be parameterized as 



where 
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/ = x/TF+Zf an d = tan/3 = /2//1. Under the SU(2)l SM gauge group, 77 is a real 
scalar, while H transforms as a doublet and can be identified as the SM Higgs doublet. The 
kinetic term in the non-linear sigma model is 
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where 5^ = ptan0 w /v 1 — tan 2 #iy/3 with 6>w being the electroweak mixing angle. As $1 
and $2 develop their VEVs, the new heavy gauge bosons Z', Y°, and X ± get their masses 
proportional to /: 
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The gauged SU(3) symmetry promotes the SM fermion doublets into SU(3) triplets. 
There are two possible gauge charge assignments for the fermions: the 'universal' embed- 
ding and the 'anomaly-free' embedding. The first choice is not favored by the electroweak 
precision data [4], so we focus on the second way of embedding. The quark Yukawa inter- 
actions for the third generation and the first two generations can be written respectively 
as 



£3 
£1,2 



i\\t c MQ 3 + iXfaQlQa + i-^d c m e ijk <$>[& 2 Q k z + h.c, 

\mn 

iXtd c ln Ql^ + i\ d 2 n d c 2n Q T n <S> 2 + i^-u c m e ijk ^* 2 j Q k n + h.c, 



i*+ c <T>t, 



(5) 
(6) 



where n — 1,2 are the first two generation indices; i, j, k = 1,2,3; Q3 = and 
Qn = {d n Li —u n L, iD nL }; d c m runs over (d c , s c , b c , D c , S c ); d\ n and d 2n are linear combinations 
of d c and D c for n — 1 and of s c and S c for n = 2; u c m runs over (u c , c c , t c , T c ). For simplicity, 
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we assume the quark flavor mixing are small and neglect the mixing effects. Eqs. and 
([6]) contain the Higgs boson interactions and the mass terms for the three generations of 
quarks: 

Ct ^ -f\\ [x{c p tl(-sit L + cxT^Gtir,) + s p t c 2 \s 2 t L + c 2 T L )G 2 (r))] + h.c, (7) 
£d n ^ -fX d 2 n [x d x n cpd c l (s 1 d nL + c 1 D nL )G* 1 (r 1 ) + s^d c 2 (-s 2 d nL + c 2 D nL )G* 2 (rj)] + h.c, (8) 
C b - -j-f 2 spcps 3 b c b L G 3 {r]) + h.c, (9) 

C g ~ -^fspCpS Z q c q L Gl{7]) + h.c. (q = u,c), (10) 



where 



51 s 8111 -^r • s ; s sm w 53 s sm v^7 ' 

GM = i + i-Lfo - !)„ - -L( t/J - 1)V, (ii) 

with h and t> being the SM-like Higgs boson field and its VEV, respectively. The mass 
eigenstates are obtained by mixing the corresponding interaction eigenstates, e.g., the mass 
eigenstates (t m i,T mL ) and (t^T^) are respectively the mixtures of (tL,T L ) and (t c ,T c ). 
The diagonalization of the mass matrix in Eqs. ([7]) and flB]) is performed numerically in our 
analysis, and the relevant couplings of h and r) bosons can also be obtained without resort 
to any expansion of v/f. Hereafter we denote the mass eigenstates without the subscript 
'm' for simplicity. 

The Yukawa and gauge interactions break the global symmetry and then provide a po- 
tential for the Higgs boson. However, the Coleman- Weinberg potential alone is not sufficient 
since the generated h mass is too heavy and the new pseudoscalar rj is massless. Therefore, 



one can introduce a tree-level \i term which can partially cancel the h mass 

- ^»t*. + h,) = -2„V%c s co S [^Jl—^ cos 
The Higgs potential becomes 



(12) 



V = -m 2 H ] H + \(H ] H) 2 - + X'H^Hrf + ■■■ , (13) 
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m = m — 



A = A 



A' 



(14) 



with mo and Ao being respectively the one-loop contributions to the h mass and the quartic 

n 

couplings from the contributions of fermion loops and gauge boson loops |4j. The Higgs 



VEV and the masses of h and 77 are given by 
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spcp \^2f S pCp 
The Coleman- Weinberg potential involves the following parameters: 



(15) 



/, x{, tp, fj,, m v , m h ,v. 



(16) 



Due to the modification of the observed W gauge boson mass, v is defined as [10] 
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where = 246 GeV is the SM Higgs VEV. Assuming that there are no large direct contri- 
butions to the potential from physics at the cutoff, we can determine other parameters in 
Eq. ( jTBl) from /, tp and m v (m^) with the definition of v in Eq. ( jTTl) . 



III. r]r], nh AND hh PRODUCTIONS AT LHC 

At the LHC the double rj production can proceed through gluon-gluon fusion and bb 
annihilation, as shown in Fig. [1] and Fig. [21 respectively. Since their Yukawa couplings are 
very small, we do not consider the contributions of qq (q = u, c, d, s) annihilation processes. 
For the gluon-gluon fusion process, there are two types of Feynman diagrams. One is the 
triangle diagrams where an off-shell (or on-shell) h boson, produced from gluon-gluon fusion 
through the heavy quark loops, decays into a pair of rj bosons. The other is the box diagrams 
where the double 77 bosons are produced through quark boxes. Due to the large Yukawa 
couplings of the heavy quarks and the large parton distribution function of gluon at the 
LHC, the contributions of the gluon-gluon fusion process can dominate over those of bb 
annihilation process. The Feynman diagrams of gg{bb) — > hh process can be obtained from 
those of gg(bb) — > rjrj by replacing the final state two rj bosons with two h bosons. 

The rjh associated production at LHC can proceed through gluon-gluon fusion and qq 
annihilation, as shown in Fig. [3] and Fig. HI In addition to the diagrams which are similar 
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FIG. 1: Feynman diagrams for double rj production via gluon-gluon fusion in the SLHM. Here 
i,j = 1,2 with (/i,/2) denoting (i,T) or (d,D) or (s,S). The diagrams by exchanging the two 
gluons or exchanging the two r/ bosons in (c,d) are not shown here. 
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FIG. 2: Feynman diagrams for double r/ production via bb annihilation in the SLHM. 

to those of 7777 production, there are some Feynman diagrams where Z and Z', produced 
from gluon-gluon fusion through the quark loops and from qq annihilation, decay into r]h as 
shown in Fig. |3]^a) and Fig. Ufa), respectively. Although the new neutral gauge boson Y° 
can also contribute to the processes gg —> r\h and qq — > rjh, its gauge coupling is suppressed 
by v/f and 1/tp @, Hi- Therefore, we neglect the contributions of Y°. 

The calculations of the loop diagrams in Fig. [1] and Fig. [3] are straightforward. Each 
.cop diagram is composed of some scalar .oop tactic Q which are ea.cu.ated by using 
LoopTools [20(. The amplitudes of triangle diagrams are as follows: 

Sab 2 / „ _ „ 

%a) + %b) = T^N^ • 4(^2) ^m/(l + 2mJC - fcj ■ k 2 C^ 1 JhJI ^ Jh '"> 
71 f 
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z h + tm h T h 



s — m 



+ 9h v ff > 



where s = A; 2 = (&i + /c 2 ) 2 . M^y) denotes the amplitude of Fig. with x = 1, 3 and 
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FIG. 3: Feynman diagrams for rjh associated production via gluon-gluon fusion in the SLHM. For 
(b-e), i,j = 1,2 with {fx, f^) denoting (t,T) or (d,D) or (s,S) ; For (a), fi denotes the SM quarks 
and new quarks T, D and S. The diagrams by exchanging the two gluons or exchanging the r] and 
h bosons in (d,e) are not shown here. 
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FIG. 4: Feynman diagrams for r\h associated production via bb annihilation in the SLHM. The 
parton q in (a) denotes u, c, d, s, b. 



y=a, b, c. Cq = Co(0, 0, s, mj, ml, mj) is the 3-point Feynman integrals scalar function, and 
a and b denote the color factor of gluons. 
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F a ^ v '(rrif, g%ff) is the effective coupling of ggZ 
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221: 



(19) 
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where g a is the coupling of axial vector current and Fi(s) ( % = 1 — 4 ) are scalar functions, 
Ft = -l(B [0,m 2 f ,m 2 f ] - B [s,m 2 f ,m 2 f ] + 1 + 2C [0,0, s,m 2 f ,m 2 f ,m 2 f ]m 2 f ), 



s 

-F 2 = F 3 = - 

s 

1 



— -(B Q [0, m 2 , m 2 ] - B Q [s, m 2 , m 2 } + 1 - 2C [0, 0, s, m 2 , m 2 , m 2 ]m 2 ) + 1 



2 

Fa = -^{B [0,m 2 f ,m 2 f ]-B [s,m 2 f ,m 2 f ] + l-2C [0,0J,m 2 f ,m 2 f ,m 2 f ]m 2 f ) + l, (22) 



where the unity in F^ is the anomaly term. 

The amplitudes expressions of box diagrams are lengthy, which are not presented here. 
The hadronic cross section at the LHC is obtained by convoluting the parton cross section 
with the parton distribution functions. In our calculations we use CTEQ6L [23] to generate 
the parton distributions with the renormalization scale /j,r and the factorization scale /ip 
chosen to be //r = \ip = 2m v for r\r\ production process (/xr = jip = m ri + m h for rjh 
production; /ir = fip = 2rrih for hh production) and use the two-loop running coupling 
constant a s with a s (mz) = 0.118. The h boson mediating gg(bb) — > r\r\ and Z' mediating 
99 (qq) ~>* vh can be respectively on-shell for > 2m ri and rriz> > m v + m^. In order to 
take into account possible resonance effects of h and Z', we have calculated the decay modes 
of h and Z', which are shown in appendix. 



The SM input parameters relevant in our study are taken as m t = 173.3 GeV 24j] and 



m z = 91.1876 GeV 



251 ] . The free SLHM parameters are /, t@, m v (mh),xf and x s x . As 



shown above, the parameters x\, /i, m h (m v ) can be determined by /, tp, m v ( mh) and v. 
To satisfy the bound of LEP2, we require that rrih is larger than 114.4 GeV [26|. Certainly, 
due to the presence of the dominant decay mode h — > r\r\ and suppression of hZZ coupling 
||], Q, the LEP2 bound on rrih should be loosened to some extent. The recent studies 
about Z leptonic decay and e + e~ — > t + t~^i process at the Z pole show that the scale f 
should be respectively larger than 5.6 TeV and 5.4 TeV, which does not depend on tp [27]. 
Such large values of / can suppress the SLHM predictions sizably. However, the factor tp 
in the couplings of h and r\ can be taken as a large value to cancel the suppression of / 
partially. For the perturbation to be valid, tp cannot be too large for fixed /. If we require 
OK/f 4 )/(D(v 2 /f 2 ) < 0.1 in the expansion of v, tp should be below 28 for / = 5.6 TeV. In 
our calculation, we take / = 5.6 TeV and tp = 15, 20, 25, respectively. 

Besides, the SLHM predicts a heavy neutrino for leptons of each generation, and the 
mixing of the heavy neutrinos with the light neutrinos in conjunction with a family mixing 
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FIG. 5: For / = 5.6 TeV, hadronic cross sections of gg(bb) —> r/r/ and gg(qq) — > rjh at the LHC 
versus the r/ boson mass. The incomplete lines for tan/3 = 20 and 15 show the lower bounds of the 
r\ mass, respectively. 



29], which can 



in the lepton sectors produces the new lepton mixing matrix |6|, 1281 
lead to lepton-flavor violating processes, such as \i — > ej, \i — > eee and /iiV — > eN. The 
experimental constraints from the three processes are very strong. For example, using only 
two lepton generations, we need / > 8 TeV or very small mixing angles or heavy neutrinos 
mass splitting, i.e., sin26 l < 0.01 or 5 < 1% with tp 



1 



However, as in the quark 



sector of the SM, this lepton-flavor violation will vanish in the limit that the mixing matrix 
Vi is diagonal or the masses of the heavy neutrinos are degenerate [6|. In this paper, we 
assume the two scenarios, so that / and tp are free from the experimental constraints of 
the lepton-flavor violating processes. Note that the parameters of lepton sector, i.e., mixing 
matrix Vi and masses of the heavy neutrinos, are not involved in our calculations directly. 

The small mass of the d(s) quark requires one of the couplings Xf and {\\ and A2) to be 
very small, so there is almost no mixing between the SM down-type quarks and their heavy 



partners. We assume Af (Af ) is small, and take x d x = 1.1 x 10 



-4 



2.1 x 10~ 3 ), which 



can make the masses of D and S be in the range of 1 TeV and 2 TeV for other parameters 
taken in our calculations. In fact, our results show that different choices of x d x and x s x can 
not have sizable effects on the result. 

In Fig. [51 we plot the hadronic cross sections of gg{bb) — > 7777 and gg(qq) — > f]h at the 
LHC versus the r\ boson mass, respectively. Fig. [5] shows the cross sections of these processes 
are all sensitive to the rj mass, and the values decrease with increasing of the 77 mass. For 
the double r] production, the cross section of gluon-gluon fusion process always dominates 
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FIG. 6: For / = 5.6 TeV, hadronic cross sections of gg — > hh and bb — > hh at the LHC versus the 
h boson mass. The pp — > hh process in SM includes gg —> hh and bb —> hh. 

over that of bb annihilation process. For tp = 25 and 35 GeV < < 75 GeV, an on-shell 
h boson can be produced from the gluon-gluon fusion and bb annihilation, and decays into 
7777. Therefore, the cross sections of the two production processes are enhanced sizably and 
reach O(10 3 fb) and O(10 2 fb), respectively. However, in the parameter space where the h 
boson is off-shell, those two cross sections can reach respectively (9(10/6) and 0(lfb) for 
tp = 25, and become smaller for tp = 20 and tp = 15. 

For the rjh associated production, the cross sections of gluon-gluon fusion process and qq 
annihilation process are about in the same order of magnitude, which can reach O(10 2 fb) 
for a light 77 boson. Our numerical results show the contributions from the process mediated 
by Z and Z' as shown in Fig. HJ^a) are dominant for qq — > rjh, but in Fig. |3]^a) are much 
smaller than other processes for gg — » t]h. Although the intermediate state Z' can induce 
the resonance effects, the contributions from the process mediated by Z can still dominate 
over those of Z' due to the large mass of Z' and the small branching ratio of Z' — > rjh. 

In Fig. [61 we plot the hadronic cross sections of gg — > hh and bb — > hh at the LHC versus 
the h boson mass. We find that the cross section of bb — > hh can be neglected compared 
with that of gg — >■ hh, and the cross sections of the two processes decrease with increasing 
of rrih- Compared with the SM prediction, the cross section of hh production at LHC in 
SLHM can be enhanced sizably for a large tp. For example, with / = 5.6 TeV and tp = 25, 
the cross section can be approximately enhanced by 80% for = 120 GeV. 

Now we briefly discuss the observable signatures of 7777, rjh and hh productions at LHC, 
respectively. For 10 GeV < < 300 GeV, the r\ boson mainly decays into bb, rf and 



10 



gg. The branching ratio of rj — > rf is about 10% of rj — > bb [10l . 1 121 ] . The huge QCD 
backgrounds make it essentially impossible to discover the signatures pp — > rjrj — > bbbb and 
pp — y TjTj — y 9999 at LHC If one of the double rj bosons decays into tt, and the other 
decays into bb, the signal to background ratio S/B can be enhanced S1Z a bly Q. Na,ne. y , 
pp — > rjrj — > rfbb is a promising channel to search for rjrj at LHC. Besides, the rare mode 
pp — > T]T] — > bb r y r ) may be also promising since the narrow 77 peak can be reconstructed 



to distinguish the signal from the backgrounds [30j. Fig. [5] shows that, for tp = 15 and 
tp = 20, the cross sections of rjrj production are only several fb and even less than 1 fb, so 
it is difficult to search for rjrj at LHC. However, for tp = 25 and 35 GeV < m v < 75 GeV, 
the cross sections can reach several thousand fb, so it is feasible to search for rjrj through 
pp — > rjrj — > rfbb at LHC. 

For rah < 150 GeV and tp = 25, the decay — >• 7777 is dominant and h — > bb is subdominant 



@, llQl, ll3|. The largest mode pp — > rjh — >■ rjiirj is interesting, and a detailed study is needed 
to establish the feasibility of searching for rjh. For rrih < 150 GeV and tp = 20 or 15, the 
dominant decay mode is h — >• bb, and h — > rjrj is forbidden kinematically. So the promising 
channel is pp — )■ rjh — >■ bbrf. The decay h — >■ W^W 7 is dominant for > 150 GeV, and the 
decay fa — > ZZ is subdominant for rrih > 160 GeV. Here the largest mode bbWW from the 
decays rj — )■ 66 and fa — )■ H^H 7 has huge background pp tt —> bbWW, so the mode is not 



optimistic 



31] . The next-largest mode bbZZ from the decays 77 — >■ 66 and fa — > may be 



more promising. 

Similar to the analysis of rjrj and rjh production processes, for rrih < 150 GeV and tp = 25, 
pp — > hh — > rjrjrjrj is the largest mode, and a detailed study is needed to analyze the signal 
and the relevant backgrounds. Besides, the promising channels are pp — > hh — > bbrf (6677) 
for m h < 150 GeV and tp = 20 or 15, and pp -> fafa ->• WWWW for 150GeV < m h < 200 
GeV. 



IV. CONCLUSION 

In the framework of the simplest little Higgs model, we perform a comprehensive study for 
the pair productions of neutral Higgs bosons at LHC, namely gg(bb) — > rjrj, gg(qq) — > rjh and 
gg{bb) — > hh. We find that the cross section of gg — > rjrj process always dominates over that 
of 66 rjrj. The cross sections of the two processes can reach respectively several thousand 



11 



fb and several hundred fb when an on-shell h boson mediates the two processes. When 
the h boson is off-shell, the cross sections of the two processes are reduced by two orders of 
magnitude, respectively. Besides, the cross sections of gg — > rjh process and qq — > rjh process 
are about in the same order of magnitude, which can reach O(10 2 fb) for a light rj boson. 
For the hh production process at LHC, the cross section in SLHM can be enhanced sizably 
compared with the SM prediction. The above results imply that it is possible to search for 
rj boson and h boson through those production processes at LHC. We briefly discuss the 
observable signatures of 7777, rjh and hh at the LHC, and a detailed analysis of the signatures 
and relevant backgrounds are necessary in the future study. 
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Appendix A: The decay widths of h and Z' 



In addition to the SM decay modes, the h boson in the SLHM can decay into tjt] and 
Zr] in the kinematically allowed parameter space. The SLHM corrections to the tree-level 
decays h — > ff, WW, ZZ are mainly from the corresponding modified couplings: 



T(h -> XX) = T(h -> XX) SM (g hXX /g™ x ) 2 , 



(Al) 



where XX denotes WW, ZZ or fermion pairs. The SM decay width F(h — > XX) sm is 
obtained using the code Hdecay [32j. ghxx and g^xx are ^ ne couplings of hXX in the 
SLHM and SM, respectively. The couplings gnww an d ghzz can be found in jlo| . 



The decay rates of h — > gg and h — > 77 are as follows 
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32n 3 v 2 



where / = t, T, D, S, Tf 

T(h -> 77) : 



4m? 



and yf 
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(A3) 
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where y v = ^rg hVV with V denoting the charged gauge bosons W ± and X ± . N c f and Qf 
are respectively the color factor and the electric charge of the fermion running in the loop. 
The dimensionless loop factors are 



F x (r) = 2 + 3r + 3r(2 - r)/(r), F 1/2 (r) = -2r[l + (1 - r)f(r)), 



(A4) 



where 



f(r) 



[sm'\l/^)] 2 , r > 1 

-j[Hv+/v-) - w] 2 , r < i 



with r]± = 1 ± yi— T. 

The decay rates of /i — >■ r/r/ and /i — >■ Z?? are 



r(/i —7- 7777) 



87r 



x 



nr. 



ml m 2 



32vr^ V " H 

where x v = 4m 2 /m 2 h and A(l, x, y) = (1 — x — y) 2 — 4xy. 

The heavy gauge boson Z' mainly decays into fermion pairs. 



2 ' 2 

ml ml 



nz' -+ fi) = ^V^r 



{i9 f L f + ( 9 f R ) 2 ) (i- T {) + \gW R T f 



m z > 



where g s L and g R are from the coupling Z' f^{g{ j P L + g R P R )f , which can be found in |6|, |33j. 



(A5) 



(A6) 



(A7) 



, — o~, — o - ; m z > 

m z , m z , 



(A8) 



where a>z> 



mz cos 9 W (1— tan 2 8 W ) 



Since the coupling is suppressed by j% and without en- 
hancement of tp, the decay mode Z' — > WW is neglected here. The decay modes Z' — > Y°Y° 



f y/3-ta.n 2 e w 
tp, the d( 

and Z' — > X + X~ are forbidden kinematically 
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